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Abstract. 

We investigate finite-size scaling of genuine multisite entanglement in the ground state of 
quantum spin-1/2 Heisenberg ladders. We obtain the ground states of odd- and even-legged 
Heisenberg ladder Hamiltonians and compute genuine multisite entanglement, the generalized 
geometric measure (GGM), which shows that for even mngs, GGM increases for odd-legged 
ladder while it decreases for even ones. Interestingly, the ground state obtained by short- 
range dimer coverings, under the resonating valence bond ansatz, encapsulates the qualitative 
features of GGM for both the ladders. We find that while the quantity converges to a single 
value for higher legged odd- and even-ladders, in the asymptotic limit of a large number of 
rungs, the finite-size scaling exponents of the same tend to diverge. The scaling exponent of 
GGM is therefore capable to distinguish the odd-even dichotomy in Heisenberg ladders, even 
when the corresponding multisite entanglements merge. 
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1. Introduction 

The sequences of both odd and even numbers mathematically reach the same infinity. The 
question is whether limits of functions of odd and even numbers reach the same limiting 
function as the respective sequences tend to infinity. In many-body physics, the dichotomy 
between the physical properties of odd and even quantum Heisenberg ladders is well known. 
For example, the seminal work in Ref. [1] shows that the spin gaps of odd and even 
Heisenberg ladders exhibit different behavior, although the values are expected to converge to 
a single limiting value, with an increase in the number of legs of the ladder. It may now 
be asked whether one can identify a physical quantity that would have different limiting 
values depending on whether odd- or even-legged ladders are followed to reach the infinite 2D 
square lattice. In this paper, we answer this question affirmatively by identifying a quantum 
information theoretic quantity that does the job. 

From the perspective of quantum many-body physics, the ground state of the Heisenberg 
ladder [1-3] is an important physical system with a rich topological order. The significance 
of these quantum spin ladders lie in their nontrivial intermediate properties between 
one-dimensional (ID) and two-dimensional (2D) spin lattices. For example, specific 
characteristics of Heisenberg ladders do not extrapolate trivially from the ID Heisenberg chain 
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to the 2D square lattice. This is due to the fact that odd and even Heisenberg ladders show 
different physical properties: Even-legged ladders are spin-gapped and have exponentially 
decaying correlation lengths while odd-legged ladders are gapless with power-law decay 
[1-5]. Investigating the odd-even dichotomy and the scaling of cooperative properties in 
large Heisenberg ladders remains an elusive proposition, primarily due to the unavailability 
of suitable analytical and numerical tools. Quantum correlations [6, 7] have been used as 
a tool to detect cooperative phenomena and topological order in ground states of Heisenberg 
ladders [8-15] . In recent years, there have been studies to understand the even-odd disparity in 
terms of entropy area law [16], Renyi entropy [17], entanglement spectra [18], etc. However 
investigating the odd-even dichotomy via the scaling of cooperative multisite properties in 
large Heisenberg ladders remains an elusive proposition, primarily due to the unavailability 
of suitable analytical and numerical tools. 

In this work, we investigate the variation in ground state properties of even versus odd 
Heisenberg ladders by analyzing its genuine multisite entanglement. To characterize the 
quantity in the ladder states, we use a computable measure, called the generalized geometric 
measure (GGM) [19] (cf. [20]). The ground state of the Heisenberg ladder Hamiltonian 
is obtained using exact diagonalization algorithms [21,22] for moderate system size. We 
observe that the genuine multisite entanglement behaves in qualitatively different ways for the 
ground states of the odd- and even-legged Heisenberg ladders - thus, detecting the odd-even 
dichotomy present in the system. In particular, the GGM increases with increasing number of 
ladder “rungs” for odd ladders while it decreases for even ladders. We subsequently observe 
that in terms of the behavior of GGM, the ground states of these ladder Hamiltonians are 
qualitatively similar to the ground states obtained from the RVB ansatz. Simulating the 
ground state of the Heisenberg ladders using RVB states allows us to analyze the finite- 
size scaling of genuine multisite entanglement in relatively large spin lattices by employing 
recursion methods [18,23,24]. We observe that although the behavior of the GGMs for 
odd and even RVB ladders are qualitatively different, they converge to a single value in 
the asymptotic limit. Therefore, for ladders with large number of rungs, as the number of 
ladder “legs” are increased, the odd versus even demarcation in terms of GGM vanishes. 
However, evaluation of the finite size scaling exponent of multisite entanglement for large 
lattices reveals that the scaling exponents tend to diverge for odd and even ladders, as the 
number of legs are increased, even though the amount of genuine multisite entanglement 
converges in the asymptotic limit. 

The paper is arranged as follows. In Sec. 2, we characterize and compute the genuine 
multisite entanglement in ground states of the Heisenberg ladder. We introduce the RVB 
ladder states in Sec. 3 and discuss the density matrix recursion method to obtain reduced 
density matrices. In Sec. 4, we compute GGM and study its finite-size scaling in RVB ladder 
states. We end with a discussion in Sec. 5. 

2. Characterization of genuine multisite entanglement in Heisenberg ladders 

The genuine multisite entanglement of a pure quantum state can be computed using the 
generalized geometric measure. For an n-party pure quantum state, \tf}n), the GGM is defined 
as the optimized fidelity distance of the given state from the set of all states that are not 
genuinely multiparty entangled: 

Gilijn)) = 1- Aliaxil^n)), ( 1 ) 

where AmaxdV'™)) = max |(x|V'ra)l with |x) being an n- party state not genuinely multiparty 
entangled. We note that a pure state is not genuinely multiparty entangled if it is product 


Diverging scaling with converging multisite entanglement 


3 


across any bipartition. For an n-party pure quantum state, |'0„), consisting of the parties Ai, 
A 2 ,. . An, Equation (1) can be shown to be equivalent to the complete form [19] 

g = l- max{A^^i|/E U L = {A,}^^^,K n L = 0}, (2) 

where Xk-.l is the maximal Schmidt coefficients in the bipartite split K : L of \tjjn)- 

We note that all possible bipartitions, K : L, of the system are considered in Eq. (2), 
with the K subsystem in the above bipartition containing all possible combinations of Ai, for 
i = 1,2,... ,n. The computation of GGM depends on the efficient generation of arbitrary 
reduced density matrices across all possible bipartition of the spin system [25]. Eor states 
where the reduced density matrices can be efficiently generated, the GGM turns out to be a 
computable measure of genuine multisite entanglement. Additional leverage is obtained if the 
state is known to be symmetric and the maximal Schmidt coefficient is known to arise from a 
selected subset of all possible bipartitions. 

Let us now consider the GGM of the ground state of spin-1/2 Heisenberg ladders, 
which have been intensively studied in strongly-correlated physics in order to investigate 
exotic quantum phenomena, like high-Tc superconductivity [3], chiral Mott insulators [26] 
etc. Such studies are also interesting in view of the fact that Heisenberg models have been 
implemented using several experimental settings, ranging from optical lattices to nuclear 
magentic ressonance [27-33]. The Hamiltonian of a quantum spin-1/2 Heisenberg model, 
with nearest-neighbor (NN) interactions, can be written as Hint = j ' 'T? ) where 

J( J >0) represents the NN antiferromagnetic (AEM) coupling constant. The indices, i and 
j, denote the sites of an arbitrary £-legged ladder, and ct] are the Pauli operators acting on 
the site. The notation \i — j \ indicates that the corresponding summation is overNN sites. 
Eigure 1, shows an £-legged ladder, with M rungs. 



Figure 1. Schematic diagram of an £-legged and Ad-rung ladder, with M (=Ad) and L (=C) 
number of spin sites along the legs and rungs, respectively. The boundary condition is shown 
by a solid line that connects the first and last sites on a specific leg. 


The model can not be analytically approached beyond ID [34]. Though various 
approximate techniques such as density matrix renormalization group [35], quantum Monte- 
Carlo [36], and RVB theory [37, 38] have been used to compute certain correlation and 
bipartite entanglement properties, the characterization of genuine multisite entanglement in 
large spin systems remains an extremely challenging task. Under these restrictions, we limit 
our exact-diagonalization study to moderate-sized Heisenberg ladders, upto 24 quantum spin- 
1/2 particles, and examine the GGM for the one-, two-, and three-legged ladders. We apply 
numerical algorithms, within the Lanczos method [21,22], in order to obtain the ground state 
of the ladder Hamiltonian, and compute the GGM. The odd- and even-legged Heisenberg 
ladders show qualitatively distinct features if one considers correlation length, energy gap 
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etc. [2]. We will now see whether such contrast in behavior can be seen by multipartite 
entanglement measure. 

2.7. Odd-legged ladders 

We now consider the GGM of the one- and three-legged (C = 1 and 3) quantum spin-1/2 
Heisenberg ladders as a function of number of rungs, A4. Fig. 2(a) displays the GGM, (/, as a 
function of number of rungs. Q exhibits alternating behavior based on whether the number of 
rungs, A4, is odd or even. This feature can not be observed in the result obtained via recursion 
technique of the RVB theory, described later in the text, as the RVB ansatz naturally requires 
an even number of rungs. From Fig. 2(a), it is evident that for both the ladders Q increases 
with increasing A4. As the number of rungs increases, the rate of increment for Q slows down 
rapidly. Interestingly note that the fluctuations of Q between odd and even rungs reduces if 
one increases number of legs which can be observed by comparing the red lines (£= 1) with 
the blue ones (£ = 3) in Fig. 2(a). 




Figure 2. Odd Heisenberg ladders: exact diagonalization vs. RVB ansatz. The behavior of 
GGM (0), with increasing even trumber of rungs (A4) for odd Zi-legged ladders, in ground 
states obtained by (a) exact diagonalization of the Heisenberg Hamiltonian and (b) using 
short-range RVB states. The solid lines show fits to the data values using Eq. (19). The 
dashed line serves as a guide to the eye. We observe similarity in behavior between the exact 
and RVB ground states with respect to behavior of Q with increasing Af. All quantities are 
dimensionless. 


2.2. Even-legged ladders 

Similarly, we also consider the GGM of the two-legged ladder {£ = 2), as a function of 
number of rungs, M. (see Fig. 3). We observe, Q decreases with the increase in M for even 
rungs, while for odd rungs, it increases. However, as seen from Fig. 3, the difference of 
GGM between even and odd rungs decreases with the increase of rungs, and for relatively 
high values of Jvi, they both correspond to a single line following the same pattern. The same 
feature is obtained for £ = \ and £ = 3 in Fig. 2(a). Therefore, we conclude that with even 
rungs, Q increases for odd legged ladders while decreases for even ones. 

At this point, a question that arises is whether the distinct qualitative features obtained 
for the GGM using exact numerical simulations of the Heisenberg model can be modeled by 
using the RVB ansatz [37,38]. This is motivated by the fact that several studies have observed 
the odd-even dichotomy in Heisenberg ladders using the RVB ground states [39]. It has 
been noticed that frustrated quantum Heisenberg spin models normally possesses short-range 
dimer states as their ground states. In particular, the ground states of the Ji — J 2 model both 
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in ID and 2D [40], the Ji — J 2 — J 3 AFM Heisenberg model [41,42], and the frustrated 
AFM on the 1/5-depleted square lattice [43], in certain parameter regimes, are the RVB 
states. Parallely, a family of rotationally invariant spin-1/2 Klein Hamiltonians exhibiting 
ground-state manifolds covered by NN valence bond states have also been proposed [44]. In 
this direction, a more systematic approach was proposed in which dimer models in different 
two-dimensional lattices like square, hexagonal, kagome, are introduced whose exact ground 
states are valence bond states [45]. Further supporting evidence for RVBs being ground states 
of Heisenberg ladders are provided in [46,47] and [4,48,49]. Recent results in the tensor- 
network formalism reveal that RVB states can be used efficiently to simulate the ground state 
properties of kagome [50] and the Ji — J 2 square Heisenberg models [51]. 



Figure 3. Even Heisenberg ladders: exact diagonalization vs. RVB ansatz. The behavior 
of GGM (Q), with increasing number of rungs (A^) for the 2-legged ladder, in ground 
states obtained by exact diagonalization of the Heisenberg Hamiltonian. The inset shows the 
behavior of Q in 2- and 4-legged ladders using short-range RVB states. The solid lines show 
fits to the data values using the equation given in Eq. (19). The dashed line serves as a guide to 
the eye. Once again, we observe the qualitative similarity between the exact and RVB ground 
states with respect to the behavior of Q with increasing even A4. All axes are dimensionless. 


In our work, we assume short-range RVB states, with NN dimer coverings, as the 
possible ground state of spin-1/2 Heisenberg ladders. Numerical investigation of the ground 
state from exact diagonalization and the RVB theory, for spin ladders upto 16 spins, provides 
considerable support for the RVB ansatz from the evaluation of the fidelity (iF) and the 
normalized relative difference in average energy (AB) [52] between the exact ground states 
and the RVB states. For example, for both the 2-leg and 3-leg quantum spin ladders, upto 16- 
spins, B as high as 0.9 and AE as low as 0.04 are obtained. These numerical findings gives us 
a good motivation for investigating the genuine multisite entanglement properties of the spin- 
1/2 Heisenberg ladder using the RVB ansatz. Let us also mention here that RVB theory has 
been also popularized as a possible theoretical tool to understand high-Tc superconductivity 
[3, 53] and are important in investigating cooperative phenomena in quantum many-body 
systems [9, 10, 54], and related to fault-tolerant quantum computation [55]. 

3. Resonating valence bond ladders 

Consider a quantum spin-1/2 ladder, with £ “legs” and j\4 “rungs” on a bipartite lattice 
(A, B), comprised of M{= A4) sites along the horizontal side and L{= £) sites the vertical 
side. The total number of spins, n (= M.N), is always even, to allow for complete 
dimer coverings. Now if the interactions between the spins are restricted to be short-ranged 
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and isotropic, we assume that only NN dimer coverings are allowed. The equal weight 
superposition of all such possible dimer coverings on the lattice would give us the so-called 
RVB state, given by 

IV'> = E [\Ai,B,) (g) |gl2, B 2 ) ® ....\An, B^)]^, (3) 

c 

where C refers to a complete dimer covering of the lattice with the summation running over 
all the coverings, and \Ai,Bj) refers to the dimer -^(| ti)| ij) — (I ii)| ti))> formed 
between spins at sites i and j, on the sublattices A and B, respectively. The RVB state \tj}) 
is rotationally invariant and is always genuinely multisite entangled state in the asymptotic 
limit [24]. 

The RVB state in Eq. (3) is unique. This is done by defining the RVB state on a bipartite 
lattice {A, B). A bipartite lattice is formulated by dividing the spin lattice into two sublattices 
A and B, such that a spin in sublattice A has spins in sublattice B as its nearest neighbours, 
and vice-versa. In our formalism, we require that all NN dimer states are directed from spins 
on sublattice A to spins on sublattice B, which removes possible ambiguity in the sign of the 
ground state, and ensures that the defined RVB state is unique. 

Let the RVB state, defined in Eq. (3), for a quantum spin ladder be denoted by |Af, £), 
with C legs and A4 rungs. Now, let us consider the system containing M = Ai + 2 spins 
along the rungs, and L = C number of spins along the legs. Eor even £, the state with open 
boundary condition can be generated recursively as [18,23,24,56] 

\M + 2,C) = |Af+ 1,£)|1) m+2 T \Ai-^ lA)\2^rn+l,m+2 
= \M^ C)\2)m+l,m+2 + \AA. — 1,£) 

^ |2)m,m+l I l)m-t-25 (4) 

where \2)m+i,m+2 and |l)m +2 correspond to the RVB ladder states, |2,£) and |1,£), 
respectively, and \2)m+i,m+2 = |2)m-i-i,m+2 - |l)m-i-i|l)m-i- 2 - Here, the subscripts denote 
the rung index. Since, for an £-legged ladder, the index £ is constant in the recursion relations, 
without loss of generality, we can remove it in the state description, so that the RVB state is 
denoted by \M). Incorporation of the periodic boundary condition leads to the following 
extension of Eq. (4) [18,24]: 

\MA2f = \M + 2)l,m-|-2 + |A^)2,m+l|2)m-|-2,l, (5) 

where all the terms on the right can be calculated by using Eq. (4) for RVB states with 
open boundary condition. Hence, and hereafter, the superscript V will indicate that periodic 
boundary condition has been used for the corresponding state. Using the recursive relation 
given in Eq. (5) we obtain the density matrix characterizing the periodic RVB ladder system, 
which is given by 

P{M+2) = PiM+2) + \M){M\2,m+l ® 12) (21 i_m-|-2 

+ (|Ad -f 2 )(Ad| 2 ,m-|-l( 2 |m-|- 2 ,l + h.c), (6) 

where, the term P(m+ 2 ) corresponds to the density matrix of the non-periodic RVB ladder, 
computed using Eq. (4). 

As mentioned earlier, our main interest lies in the multisite entanglement properties 
of these RVB ladders. In order to explore this, we first need to have expressions for all 
possible reduced density matrices of the system. The maximal Schmidt coefficients obtained 
from these reduced density matrices allow us to compute the GGM of the RVB ladder, using 
Eq. (2). As the number of spins in the RVB ladder increases, there is a rapid growth of the 
number of possible reduced density matrices. However, the symmetry of a periodic RVB state 
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can be exploited to obtain the maximal Schmidt coefficient, which is required to compute 
GGM in Eq. (2), without considering all possible reduced states. For example, extensive 
numerical studies upto 16 spins confirm that for an \M^C) ladder, optimization over the 
restricted set of all reduced density matrices contained within a reduced 2 x £ block, say at 
sites m +1 and to + 2, is sufficient to obtain the maximum Schmidt coefficient for calculating 
the GGM. The symmetry and periodicity of the RVB ladder ensures that all reduced 2 x £ 
block, between any adjacent pair of sites in the lattice, are topologically equivalent. This 
reduces the computational difficulty in calculating the genuine multisite entanglement as the 
optimization over all reduced states is now limited to a 2 x £ block, which can be analytically 
derived using a recursion method as discussed in the following segments. 

For an RVB ladder with open boundary, the reduced density matrix of a 2 x £ block is 
obtained by tracing out all the spins except those at the rungs to + 1 and to + 2, as given by 

P{m+l,m+2) =-Vm |2)(2|(m+l.m+2) T-Vm-IPto+I 

® |l)(l|(m+2) + (|2)m+l,m+2(l|m+2(XAt|m+l + h.C.) (7) 

where AOvi = {M\M) and 

Pm+l = tfm (|2) (2|77i,m+l )? und (8) 

{XM\m+l = — 1\M). (9) 


By using Eq. (7) we obtain the reduced density matrix for the 2 x £ block at rungs m +1 
and TO + 2, for the periodic RVB ladder state, P^_|_ 2 )> given by Eq. (6). The reduced density 
matrix is given by 

P(rn+l,m+2) ~ Pm+l,m+2 + b.m+l 12) (21 1 

+ {\J^) 2,m+l\^) l,m+2\{J^ + 2| + h.C.), 

= Pm+l,m+2 + ^m+l,m+2 + (Cm+l,m+2 + h-C-)i (10) 

where 


?m+l,m+2 — A/Od —l|l)(l|m+l ® Pm+2 + ■N'M — 2Pm+l 
® Pm+2 + (|l)(XA^-l|m+l) <8) Pm+2 + h.C.), 




2 _ 
m+l,m+2 


X 

X 

X 


|2)m+l,m+2 {1 |m+l (XAI 1171+2 
M 

i |m+2 {XjM —i |m+l + Pm+l ® |l)m+2 

1 

{XM\m+2 + l/A/l(|/Ci)m+l|l)m+2) 

M 

(l|m+l'y^(^|m+2+V(-i- 

i=l 


Here — (Ad|Ad 1) and (/Ci|m+i — m,m+i(2|/C2_i)m with |/Co)m 
recursion relation for the inner product {M\M.) can now be expressed as 


( 11 ) 


( 12 ) 
|l)m. The 


A/+= A/ij\/+_i + A/’2 A/a/(_2 + 2(—1)"* (13) 

i 

where A/^ = {2|2) and all the 7 i’s can be calculated using the linear equation 

{7*|2)m,m+l = (-1)™“^^^ . 7il7i)m+l, 


(14) 
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where | 7 _, ) form an independent set of vectors consisting of certain singlet combinations of 
an (1, A/" + 2 ) spin system, e.g I 71 ) = |1). 

For odd C, the recursion relations are much simpler as the number of possible coverings 
is lower. The recursion relation for the RVB ladder with periodic boundary conditions is given 
by 

\A4 + 2, £)'P = |Af, £)l,m|2)m+l,m+2 

+ |Ad, £)2,m+l|2)m+2,l- (15) 

The reduced density matrix for the 2 x C block, at sites m + 1 and m + 2, corresponding to 
the above state is given by 

P{m+l,m+2) ~ -^M |2) (2|m+l,m+2 + J^M — 2Pm+l ® Pm+2 

+ (|2)m+l,m+2('C^|m+l,m+2 + h.C.) , (16) 

where 

|m+l,m +2 (211,771+2 {M\ 2 ,m+l and (17) 

Pm-\-l — trjTj (|2) ) (18) 

Hence, using Eqs. (10) and (16) for even- and odd-legged RVB ladders, respectively, one 
can obtain the reduced density matrices for the 2 x C block necessary to compute the 
generalized geometric measure. We note that the maximal Schmidt coefficient is obtained 
by considering the reduced states within the 2 x C block. It is observed that the maximum 
Schmidt coefficients are typically obtained from the 2 : rest or the 4 : rest bipartitions where 
the reduced spins are nearest neighbors, though there does not seem to be any distinctive 
pattern that can systematically differentiate between the typical bipartitions in odd and even 
ladders. Moreover, no systematic pattern is observed in which topologically inequivalent 
reduced states provide the same maximum Schmidt coefficient. 

We subsequently compare the GGM of the RVB state with that of the ground state of 
the Heisenberg ladder obtained by exact diagonalization. Note here that although we use 
the above method for calculating the GGM, the same recursion can be used to calculate 
other system properties like magnetization, susceptibility, classical correlators, bipartite 
entanglement and other quantum correlations, etc. A more developed exposition and 
formalism for the density matrix recursion method can be obtained in Refs. [18,24]. 

4. Diverging scaling with converging multisite entanglement 

Applying the recursion technique, we can investigate the behavior of genuine multisite 
entanglement of the RVB state in large quantum spin lattices. For example, one can study 
the finite-size scaling of GGM in £-legged ladders with large number of rungs and investigate 
the odd-even dichotomy in the asymptotic limit. For odd-legged RVB ladders, the GGM 
initially increases with increasing number of even rungs, Ai, before approaching a constant 
value at large Ai (see Fig. 2(b)), while for even-legged ladders, the GGM decreases with 
increasing A4, before flattening to a constant for larger number of rungs as shown in the inset 
of Fig. 3. Note that the DMRM approach is not possible to access an odd number of rungs. 
Importantly, we find that the behavior of genuine multisite entanglement of the ground state 
of the Heisenberg ladder with even- and odd-legged ladders is qualitatively similar to results 
obtained with the RVB ansatz. This is clearly seen by comparing Fig. 2(a) with 2(b), and the 
main with inset in Fig. 3. 

The similarity between the two methods, viz. exact diagonalization of the Heisenberg 
ladder and RVB ansatz on the same lattice, motivates us to perform finite-size scaling analysis 
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of GGM, by using RVB theory, wherein we can handle large lattice sizes. The analysis would 
shed light on the finite-size behavior of multiparty entanglement of the original Heisenberg 
ladder. The finite-size scaling of GGM in a pure quantum spin ground state, \ijj), can be 
analyzed through the scaling relation, GHip)) ~ Gci\'>P))ikn~^, where, n is the total number 
of spins, Gc is an estimated value of GGM at high n, and x is the “scaling exponent” with 
which the GGM approaches its asymptote at large n. For an £-legged RVB ladder, written as 
|£, M.), the finite-size scaling is given by the relation, 

G{\C,M))^GciC)±kn-^^^\ (19) 

Using the DMRM method, we have computed the GGM for RVB ladders upto £ < 8, with 
Ai = 20. Once can easily extend the computation for higher values of Ad. 



Figure 4. Diverging scaling with converging multisite entanglement. The behavior of the 
asymptotic GGM {Qc{C)) and the scaling exponent (x{C)) with increasing C,. We obseiwe 
that even though, the Qci^C) for odd- and even-legged ladders, converge with increasing C, the 
scaling x{C) diverge. All axes are dimensionless. 


Fig. 4 shows the values of Gc{C.) and x{C), for different values of £, where the GGM 
is scaled upto Ad = 20 rungs for an RVB ladder of C legs. We observe that as C increases, 
the Gc{C) for the odd- ladders converges to that for the even ones. This is consistent since 
the pseudo-2D spin ladders slowly approaches the square-2D lattice, and in the asymptotic 
limit, one cannot distinguish whether the system was originally generated by increasing C in 
an odd- or even-legged ladder. However, we find that the scaling exponent, x{C), for odd 
and even ladders, converges to different values with increase of C (see inset of Fig. 4). We 
therefore have a diverging scaling exponent for odd- and even-legged ladders, even though 
the corresponding multisite entanglement converge. The diverging x{C), therefore, shows 
that the finite-size scaling of GGM for RVB ladders can highlight the odd-even dichotomy at 
large C. The results show that the GGM for odd RVB ladders converges slower than that for 
even ladders at low £, which is reversed as C is increased and holds even at large £, where 
Gc{C) for odd and even RVB ladders are indistinguishable. 

Therefore, we observe that although the value of genuine multisite entanglement can 
not distinguish the odd-legged ladders from the even-legged ones for large lattice size, the 
corresponding finite-size scaling exponents are capable of detecting the difference. 

5. Discussion 

To summarize, we investigate the behavior of genuine multipartite entanglement of the 
ground state in odd- and even- Heisenberg ladders. Even though such models have immense 
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fundamental and practical importance, owing in particular to the dissimilarities on the two 
sides of the odd-versus-even divide, they are not analytically accessible. In this work, we 
began our investigation through exact diagonalization techniques to find that the genuine 
multisite entanglement, as quantified by the GGM, of the ground state obtained from the 
odd-legged ladder, increases with the number of rungs. The opposite is true in the even¬ 
legged ladder. This feature is in good qualitative agreement with the assumption that ground 
states of odd and even Heisenberg ladders are RVB states. We perform scaling analyses of the 
RVB states on ladders of large system sizes by employing the DMRM, and find that while the 
GGM of the RVB states on large ladders converges to a single value independent of the odd- 
even parity of the ladders, their scaling exponents diverge from each other. While the study 
reported is for the isotropic Heisenberg model, we have carried out parallel studies for the 
quantum XXZ model. We observed qualitative similarity of the results obtained for values 
of the zz vs. XX anisotropy up to approximately 1.4. 
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